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Abstract 

A search for supersymmetry involving events with at least one photon, one electron 
or muon, and large missing transverse momentum has been performed by the CMS 
experiment. The data sample corresponds to an integrated luminosity of 19.7 fb 1 of 
pp collisions at i/s = 8 TeV, produced at the CERN LHC. No excess of events is ob¬ 
served beyond expectations from standard model processes. The result of the search 
is interpreted in the context of a general model of gauge-mediated supersymmetry 
breaking, where the charged and neutral winos are the next-to-lightest supersymmet¬ 
ric particles. Within this model, winos with a mass up to 360 GeV are excluded at the 
95% confidence level. Two simplified models inspired by gauge-mediated supersym¬ 
metry breaking are also examined, and used to derive upper limits on the production 
cross sections of specific supersymmetric processes. 
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1 Introduction 

The extension of the standard model (SM) of particle physics through the concept of super- 
symmetry (SUSY) 111, which imposes a symmetry between fermions and bosons, can offer a 
solution to some of the issues not accommodated in the SM, such as the existence of dark mat¬ 
ter in the universe or the extreme fine tuning required to control radiative corrections to the 
Higgs boson mass (hierarchy problem) [2]-Sl- The minimal supersymmetric standard model 
(MSSM) EHZl provides a calculable framework with a fully known particle content, introduc¬ 
ing a superpartner for each SM particle. For example, squarks, gluinos, and gravitinos are the 
SUSY partners of quarks, gluons, and gravitons, respectively The MSSM has guided the search 
program for physics beyond the SM at facilities such as the Fermilab Tevatron and CERN LHC. 
Existing searches have not yet found evidence for SUSY, but a large parameter space of the 
MSSM remains to be explored. 

Within the MSSM, scenarios based on gauge-mediated SUSY breaking (GMSB) 118141811 are of 
particular interest because of their ability to naturally circumvent the so-called SUSY flavour 
problem H19I . The framework of general gauge mediation (GGM) If20l offers a clear definition 
of GMSB and establishes its key aspects. For example, GMSB predicts the gravitino (G) to be 
the lightest supersymmetric particle (LSP). The combination of this feature and the weakness of 
the coupling of G to other MSSM particles has specific consequences in collider phenomenol¬ 
ogy. Under the assumption that R-parity 10 is conserved, SUSY particles are pair-produced at 
the LHC. Except for direct LSP pair production, each SUSY particle initiates a decay chain that 
yields the next-to-lightest supersymmetric particle (NLSP). Branching fraction for the SUSY 
particle decay involving G is negligible except for the NLSP, leaving the decay of the NLSP to 
its SM partner and the G as effectively the only gravitino production mechanism. The grav¬ 
itino escapes detection, leading to missing momentum in the event. The signature of a GMSB 
signal is thus strongly dependent on the identity of the NLSP. In most GMSB models, the NLSP 
is taken to be a bino- or wino-like lightest neutralino, where a bino and wino are the super¬ 
partners of the SM U(l) and SU(2) gauge fields, respectively. Previous searches for a GMSB 
signal typically exploited the diphoton signature If2ll429l| . in which each of the two bino-like 
neutralinos decays promptly into a photon and a gravitino. Similar scenarios with nonprompt 
NLSP decays have also been considered Il30l (3l!l . Thus far, no evidence for GMSB SUSY has 
been observed, resulting in upper limits on the production cross sections given as a function of 
the SUSY particle masses, the NLSP lifetime, or other model parameters. 

This paper presents a search for SUSY with the CMS experiment at the LHC, and targets GGM 
models with wino-like NLSPs. The data sample corresponds to an integrated luminosity of 
19.7 flU 1 of pp collision data collected in 2012 at y/s = 8 TeV. In particular, we study the 
wino co-NLSP model H32l , in which nearly mass-degenerate charged and neutral winos are 
significantly lighter than the other electroweakinos and constitute the lightest SUSY particles 
aside from the gravitino. Although the lifetime of the NLSP is effectively a free parameter in 
GGM phenomenology, a prompt decay of winos is assumed in this analysis. A signature of at 
least one photon ( 7 ), one electron or muon (£), and large missing transverse momentum (pj' ss ) 
is used in this search. The photon is assumed to be emitted by a neutralino NLSP, and the 
leptons by either a charged or neutral NLSP decaying to a W or Z boson, respectively. This 
signature suppresses many SM backgrounds, obviating the need for additional requirements 
such as associated jet activity. The diagrams in Fig. [T| provide examples of the decay chains 
studied in this analysis. The present search is sensitive to the direct electroweakino production 
mode of Fig. [l](a), where the winos are produced without involving coloured SUSY particles, 
but also to strong production modes such as the gluino (g) pair-production process shown in 
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2 CMS detector 


Fig. [l](b). Similar searches were conducted by the ATLAS f33l and CMS [34:, [35] experiments 
using LHC pp collision data at = 7 or 8 TeV, as well as the CDF experiment If36l at the 
Tevatron using pp collision data at i/s = 1.8 TeV. None of these analyses sees an excess of 
events over the respective SM predictions. The wino co-NLSP model has also been probed 
through the signatures of three leptons or two leptons and two jets H37l [38 11, which target the 
decay of the neutralino NLSP to a gravitino and a Z boson rather than to a gravitino and a 
photon. None of these analyses observed a significant excess of events over their respective 
SM predictions. 



Figure 1: Diagrams showing the production and decays of wino-like co-NLSPs (y, and yj) 
leading to final states with a photon, an electron or muon, and missing momentum from unde¬ 
tected gravitinos G (a) without and (b) with involvement of coloured SUSY particles. 


2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintilla¬ 
tor hadron calorimeter (HCAL), each consisting of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. A detailed description of the CMS detector, together with a definition of 
the coordinate system and the relevant kinematic variables, can be found in Ref. [|39| . 

In the barrel section of the ECAL, an energy resolution of about 1% is achieved for unconverted 
and late-converting photons with transverse energy Ej ~ 10 GeV. The remaining barrel pho¬ 
tons have a resolution of about 1.3% up to a pseudorapidity \rj\ < 1.0, rising to about 2.5% for 

M = 1-4 BQl. 

The electron momentum is determined by combining the energy measurement in the ECAL 
with the momentum measurement in the tracker. The momentum resolution for electrons 
with transverse momentum pj ~ 45 GeV from Z —> e 1 e decays ranges from 1.7% for non¬ 
showering electrons in the barrel region to 4.5% for showering electrons in the endcaps Il4l1l . 

Muons are measured in the range \t]\ < 2.4, with detector elements based on three technolo¬ 
gies: drift tubes, cathode strip chambers, and resistive plate chambers. Through the matching 
of track segments measured in the muon detectors with tracks measured in the tracker, a trans¬ 
verse momentum resolution of 1.3-2.0% is achieved for barrel muons with 20 < pj < 100 GeV. 
In the endcaps, the resolution increases up to around 6%. The pj resolution in the barrel is 
better than 10% for muons with transverse momentum up to 1 TeV C42H . 
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Physics objects are defined using the particle-flow (PF) algorithm 11431144 1. which reconstructs 
and identifies individual particles through an optimized combination of information from dif¬ 
ferent elements of the CMS detector. The PF candidates are classified as photons, charged 
hadrons, neutral hadrons, electrons, or muons. Finally, the CMS detector is nearly hermetic, 
permitting accurate measurements of Pj' ss . 


3 Data collection and event selection 

The search is conducted in the electron-photon (e 7 ) and muon-photon (py) channels. The data 
samples are collected using a dedicated trigger for each channel, as described below. An event 
is considered to be in the e 7 (py) channel if it contains at least one high-energy photon and 
an electron (muon). Events that simultaneously satisfy the criteria for the two search channels, 
representing about 0.1% of the selected events, are classified as p 7 candidates because muon 
objects are less often the result of hadron misidentification than are electron objects. 

The trigger for the e 7 channel requires at least two isolated photon-like objects, with Ej thresh¬ 
olds of 36 and 22 GeV for the highest and second-highest Ej photon, respectively The trigger 
does not veto photon objects that can be matched to a track, allowing events with a photon 
and an electron to also satisfy the trigger. The pp channel uses a muon-photon trigger with a 
px threshold of 22 GeV for both the photon and muon objects. To ensure a fully efficient trig¬ 
ger and a similar selection efficiency for the two channels, the subsequent analysis requires 
Ex > 40 GeV for the photon and px > 25 GeV for the electron or muon. With these require¬ 
ments, the trigger efficiency for the signal models described in Section [7] is found to be 93-98% 
for both channels, depending on the model and SUSY mass values. 

Photon candidates are reconstructed from clusters of energy in the ECAL Ii40l . The momentum 
vector of the photon points from the primary pp interaction vertex to the center of the ECAL en¬ 
ergy cluster, under the assumption that the photon originates from the primary vertex, which 
is defined as the vertex with the highest X] Pj of associated tracks. Only photons from clus¬ 
ters in the pseudorapidity range \t]\ < 1.44 are included in this analysis. These clusters were 
selected as photon candidates by a set of criteria that are designed to achieve a 90% identifi¬ 
cation efficiency for true photons. For a cluster to be identified as a photon, its shape must 
be consistent with that expected from a photon, and the energy detected in the HCAL behind 
the cluster cannot exceed 5% of the ECAL energy To further suppress the misidentification of 
hadrons as photons, a PF-based isolation requirement is imposed. The transverse component 
of the momentum sum of each of the PF photons, charged hadrons, and neutral hadrons within 
a cone of A R = \J (A rj) 2 + (A <p) 2 = 0.3 around the direction of the photon candidate (where <p 
is the azimuth measured in radians) is required not to exceed fixed values defined to achieve a 
desirable balance between the identification efficiency and misidentification rate. The photon 
object that is being identified is not included in the isolation sums, and charged hadrons are 
included only if they are associated with the primary vertex. The px sums are corrected for 
contributions from additional pp interactions (pileup). To distinguish photon candidates from 
isolated electrons, photon objects are vetoed if a matching track segment from the inner tracker 
is identified. 

Electron (muon) candidates must lie in the pseudorapidity range \t]\ < 2.5 (2.4). For electrons, 
the transition region 1.44 < \rj\ < 1.56 between the barrel and the endcap detectors is vetoed 
because the reconstruction efficiency in this region is difficult to model. Electron objects are 
reconstructed by associating a cluster of energy deposited in the ECAL with a reconstructed 
track. The electron selection B45l is based on the shower shape, the matching of a track to 
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3 Data collection and event selection 


Table 1: Summary of event selection requirements and observed number of events after apply¬ 
ing the listed selection requirements in successive order. The symbol rn e - r and m z denote the 
invariant mass of the electron-photon system and the nominal Z boson mass, respectively 


Selection requirement 

e 7 channel 

pj channel 

Trigger 

26 733 051 

19456 571 

>1 accepted j 

2 718 364 

243 664 

>1 accepted l {£ = e, p) 

70 736 

32173 

AR(j,£) > 0.8 

68168 

30 232 

\m ei — m z \ > 10GeV 

29169 

— 

Ef ss > 120 GeV, M t > 100 GeV 

110 

152 


the cluster, and isolation, where the isolation variable is calculated from the momenta of PF 
photons, charged hadrons, and neutral hadrons within a cone of A R = 0.3 around the electron 
direction, corrected for the effects of pileup. The isolation sum is required not to exceed a fixed 
fraction of the electron pj, where the selection criteria are defined to obtain an 80% electron 
identification efficiency. The muon selection Il42l . targeting a 90% efficiency for true muons, 
utilizes the quality of the track fit, the number of detector hits used in the tracking, and the 
isolation. The isolation requirement for muons is similar to that for electrons, but uses a larger 
cone size A R = 0.4. Electrons and muons must originate from a primary vertex, with respective 
distances of closest approach for electrons (muons) of less than 0.2 mm (2 mm) in the transverse 
plane and <1 mm (<5mm) along the beam direction. 

The reconstruction of jets and Py 1iss is also based on the PF objects. All reconstructed PF can¬ 
didates are clustered into jets using the anti-fcj clustering algorithm [j46i [471, with a distance 
parameter of 0.5. Jet objects are used to calculate the Hj variable, defined as the scalar pj sum 
of jets. To be considered in the Hj sum, a jet must have a calibrated and pileup-corrected Il48l 
pj value greater than 30GeV, \t] < 2.5, and be consistent with an origin at the primary ver¬ 
tex i49il . In addition, it must be no closer than A R = 0.5 to the photon or lepton candidates. 
The missing momentum p| 11iss is given by the negative of the vector pj sum of all PF objects, 
with jet-energy corrections applied. The magnitude of p™ iss is referred to as Ey ISS . 

To suppress the background from final-state radiation events with an on-shell W (Z) boson that 
decays to ivj (ttj), the highest-Ej photon in an event must be separated by A R > 0.8 from 
the highest-px electron or muon. Additionally, for the e 7 channel, the invariant mass of the 
electron-photon system is required to differ by more than 10 GeV from the nominal Z boson 
mass Il50l . to reduce background from electrons misidentified as photons. 

After applying the selection requirements described above, the obtained event yields are com¬ 
pared to expectations from SM background processes. The signal region of interest is defined 
by Ey ISS > 120 GeV and Mj > 100 GeV, where transverse mass Mj is defined by Mj = 
VlEf ss p c T [l — cos A<p(£, p™ iss )], with p y the transverse momentum of the highest-px lepton 
and A(p(£,p™ iss ) the azimuthal angle between the lepton and p™ iss . The Mj requirement re¬ 
duces backgrounds from processes that produce W bosons. Table [l] shows the observed num¬ 
ber of events at different stages of the selection process. Because of a higher selection efficiency 
for muons, after implementing the selection requirements, the number of observed events in 
the pj channel is larger than in the e 7 channel. 
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4 Background estimation 

Three sources of SM background are considered: misidentified photons, misidentified leptons, 
and electroweak backgrounds. 

4.1 Misidentified-photon background 

The background from misidentified photons arises from events in which a photon object does 
not correspond to a genuine prompt photon. The dominant background processes in this cat¬ 
egory are Drell-Yan dielectron (qq —> 7 * —> e ! e ) and W (—? /Yj+jets production, in which 
an electron or jet, respectively, is misidentified as a photon. Minor contributions arise from tt 
events with leptonic top quark decays, for both the ej and f /7 channels. Events with tt produc¬ 
tion also contribute to the background if a jet is misidentified as a photon. An electron can be 
misidentified as a photon if it fails to register track seeds due to detector inefficiencies such as 
non-operational sensors in the tracker. A jet can be misidentified as a photon if a large fraction 
of its energy is carried by mesons decaying to photons, such as n° —> 77 . These two types of 
background are estimated from data using weighted control samples. The method proceeds in 
two steps. First, a control sample enriched in particles that are prone to be misidentified as pho¬ 
ton candidates, i.e., electrons or neutral hadrons, is selected by inverting certain criteria in the 
photon identification, while keeping the other selection requirements identical to those for sig¬ 
nal candidates. This control sample is called the proxy sample. The second step is to determine 
the transfer factor N misld /N proxy , where A/ mis,d is the estimate of the number of misidentified 
events in the signal candidate sample and N proxy is the number of events in the proxy sample. 
The proxy sample is then scaled by the transfer factor. The definition of the proxy sample is 
tuned to make its kinematic properties similar to those of events with misidentified photons. 
Thus, this two-step procedure takes the set of misidentified events in a control region where 
the SUSY signal contribution is expected to be negligible, e.g. in events with small E^ 1ISS , and 
utilizes it to model the distribution of the misidentified background for a given kinematic vari¬ 
able. In particular, from the extrapolation of the observed events in the control region, the 
method predicts an expectation for the number of events and corresponding kinematic distri¬ 
bution in the signal region. A detailed validation of this background estimation is performed 
by applying the method to Monte Carlo (MC) simulation samples, and comparing the outcome 
of this procedure to the known generated MC content. Good agreement is found in all such 
tests. 

The proxy sample for events with an electron-to-photon misidentification is constructed by 
inverting the electron-seed veto in the photon identification, which turns the photon candi¬ 
date into an electron proxy. The transfer factor for this proxy sample is determined by count¬ 
ing Z —» e ' e decays in a separate control sample, defined by E™ ISS < 70GeV. The ratio 
jymisid y jyP rox y cons titutes the transfer factor, where Nj) lisid is the number of Z —> e 1 e events 
in the control sample with an e + or e misidentified as a photon, while Ny U ' xy is the number of 
Z —t e + e~ events in the control sample with the proxy condition applied. The control sample 
for the e 7 channel is taken from the data set collected with the same diphoton trigger as the 
signal candidates, while the sample for the ^7 channel is from a data set based on a trigger 
that requires at least one isolated electron. To ensure that the samples dominantly consist of 
Z —» e 1 e decays, events with one high-purity electron object ("tag") are selected from the 
respective data sets. The photon candidate and the electron-proxy object are called "probes". 
For each sample of probe candidates, a fit is performed to the invariant mass distribution of the 
tag-probe system to extract the number of Z —> e 1 e decays. 

The "tag-and-probe" method described above H5ll is executed in bins of three variables: the 
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4 Background estimation 


transverse momentum of the probe object {pj Tobe ), the track multiplicity of the primary vertex 
(N trac k)/ and the number of reconstructed interaction vertices in the event (N vtx ). To account for 
the correlations in the distributions of the three variables, the dependence of the transfer factor 
on these quantities is modelled by a three-dimensional parametric function, which is then used 
to assign an event-by-event weight to the proxy sample. The transfer factor is a decreasing 
function of p^ lobe and N tiac ^, and an increasing function of Nvtx- For a median value of N vtx , 

its value varies from 0.04 for events with low p^ robe and low N trac k, to 0.007 for high p^ lobe 
and high N trac ^. The relative uncertainty in the transfer factor is typically of order 10%, which 
is dominated by systematic uncertainties such as those arising from the tag-and-probe fitting 
procedure and the parametrization of the transfer factor. The dependence of the transfer factor 
on Nvtx is approximately linear, with a value that changes from about 0.005 to 0.012 for N v t x 
values between 5 (low pileup) and 25 (high pileup). 

The estimation of the jet-to-photon misidentification background follows the same procedure 
of defining a proxy sample and scaling it with the transfer factor. The proxy sample for events 
with a jet-to-photon misidentification is constructed by inverting the requirements on the vari¬ 
able describing the ECAL cluster shape (u w in Ref. EU) and on one of the isolation variables in 
the photon selection. The transfer factor for the hadronic-proxy sample is determined through 
an assessment of the fraction of events with jet-to-photon misidentification among the photon 
candidates. This fraction is denoted as the "hadron fraction". This measurement is performed 
in a low-£" llss control sample from a fit to the distribution of a ^ based on two templates, one 
representing pure photons and one modelling the events with jet-to-photon misidentification. 
The fit is performed with photon candidates in muon-photon events, where a very small con¬ 
tamination of misidentified electrons is expected in the photon sample. The pure photon tem¬ 
plate is obtained from Z —> puy data by tagging two muons and requiring the three-body lift'} 
invariant mass to be consistent with the Z boson mass. The template that models events with 
jet-to-photon misidentification is obtained by inverting the isolation requirement on the signal- 
photon candidates. The hadron fraction is measured in pj bins of the photon candidate and 
decreases, in general, as a function of pj. In the ey channel, its value varies from 0.25 ± 0.03 at 
p x = 40 GeV to 0.08 ± 0.02 at pj = 120 GeV. In the py channel, the corresponding values are 
0.30 ± 0.03 and 0.09 ± 0.02. The uncertainties are dominantly due to possible mismodelling of 
the fit templates. The small difference in the hadron fraction of the photon candidates between 
the ey and py channels originates from small differences in trigger requirements on the photon 
object between the diphoton and muon-photon triggers used to select the ey and py data sets. 

The pj distribution of the photon objects is multiplied by the hadron fractions determined as 
described above. In the ey channel, the estimated pj distribution of misidentified electrons is 
subtracted first. The resulting distribution provides the estimate of the pj shape for the events 
with jet-to-photon misidentification. Rather than forming the ratio of this distribution with 
the pj distribution of the proxy sample, both distributions are parameterized individually by 
simple analytic functions. The ratio of these two parameterizations constitutes the transfer 
factor for the jet-to-photon misidentification. 

4.2 Misidentified-lepton and electroweak backgrounds 

The misidentified-lepton and electroweak (EWK) backgrounds are evaluated together, as de¬ 
scribed below. A misidentified lepton is defined as a reconstructed lepton that does not arise 
directly from W or Z boson decays, nor from r decays that originate from a W or Z boson. 
Misidentified-lepton events arise primarily from the decay of heavy-flavour quarks and from 
hadrons misidentified as leptons, with other sources such as decays-in-flight constituting a 
much smaller contribution. Events where both the lepton and photon are misidentified, which 
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constitute up to 30% of the total misidentified-photon background, are already accounted for 
by the procedure described in Section 4.1 The SM electroweak background is dominated by 


events with Vy (V=W, Z) production. In particular, W 7 events have the same signature as sig¬ 
nal events: an energetic photon, a lepton, and significant E“ 1SS . The EWK background includes 
rare multiboson events and events with tt 7 production, which we collectively refer to as the 
"rare EWK" background. Rare EWK events provide only a minor contribution to the overall 
background but are relevant in the high-E™ ISS signal region. 


Similar to the determination of the misidentified-photon background, proxy samples are formed 
and scaled by transfer factors to estimate the contribution of misidentified leptons to the signal 
region. Each event in the proxy sample contains at least one candidate photon and at least one 
misidentified-lepton proxy, but no candidate lepton. Proxy objects that model misidentified 
leptons are selected by inverting the isolation condition in the lepton selection. For electrons, 
the track-cluster matching requirements are also inverted to further enrich the proxy sample 
with hadronic objects. The calculation of the transfer factor used to evaluate the misidentified- 
lepton background is described below. 


The modelling of the EWK background is based on MC simulation. Samples of W 7 , Z 7 , tt 
7 , and WW 7 events, listed in the order of decreasing overall background contributions, are 
generated with up to two additional partons using the MADGRAPH5 1.3 H52| event generator. 
The PYTHIA 6.4 H531 program is used to describe the parton shower and hadronization. The 
pythia program is further used to generate samples of WZ 7 events, which produce an even 
smaller background contribution than WW 7 events. All samples use the CTEQ6L1 Il54l parton 
distribution functions (PDF). Simulated minimum-bias events are overlaid on the main hard- 
scattering events to simulate pileup. The generated particles are processed through the full 
CMS detector simulation framework based on the Geant4 Il55l package, and are subjected to 
the same event selection procedure as the data, including the trigger requirements. Differences 
between simulation and data in the pileup profile, trigger efficiency, and object identification 
efficiency are corrected by reweighting the MC events by factors that lie within a few percent 
of unity. The tty, WW 7 , and WZ 7 samples are normalized to the integrated luminosity of the 
data using cross sections calculated with the event generators, which are valid to leading order 
(LO) in quantum chromodynamics. For the tty sample, a next-to-leading order (NLO)-to-LO 
scale factor of 2.0 l56| is applied to the cross section to account for higher-order contributions. 


For the Vy background, calculated cross sections are used to fix the ratio between the Wy and 
Zy components, but the overall normalization of the combined sample is derived from data 
to mitigate potential uncertainties in the theoretical calculation. This is accomplished through 
a two-component template fit describing the Vy and misidentified-lepton backgrounds. The 
templates originate from two background samples obtained using the event selection criteria 
for the Vy MC sample and for the misidentified-lepton proxy sample. Distributions of the 
variable A (p(l,p™ ss ) from the two background samples in the control region 40 < Ej ,ss < 
70 GeV are employed as templates. The lower bound Ej ,ss = 40 GeV is applied to reduce the 
contribution of Zy events. Expected contributions from the misidentified-photon background 
and rare EWK backgrounds are subtracted from the data distribution before the fit. The fit 
provides scale factors for the template histograms that are used directly as transfer factors 
for the Vy and the misidentified-lepton proxy samples. Besides avoiding a reliance on the 
value of the theoretical Vy cross section, which is observed to underestimate the measured 
production rate of Wy events H57II58T this method has the benefit that it does not double count 
the contributions of background events with both a misidentified photon and lepton. This 
class of events is already accounted for in the misidentified-photon background sample, as 
mentioned above. 
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5 Systematic uncertainties 


Figure [2] shows the results of the template fit, which is performed in the e 7 and p r y channels 
independently. The resulting scale factors for the V 7 background in the e 7 and f /7 channels 
are 1.59 ± 0.27 and 1.47 ± 0.16, respectively, which are similar to each other as expected. The 
uncertainties in the scale factors are estimated through toy MC studies repeating the fit after 
changing the contributions of the subtracted misidentified-photon and rare EWK components 
by their uncertainties. 


19.7 ft)' 1 (8 TeV) 


19.7 fb' 1 (8 TeV) 




Figure 2: Results of the A <p(£, pj' ss ) template fit for events with 40 < £™ ISS < 70GeV, used to 
determine the V 7 and misidentified-lepton background for the (a) e 7 and (b) 1/7 channels. 


5 Systematic uncertainties 

Table [2] summarizes the sources of systematic uncertainty for the background predictions and 
the signal yields. For each source, the size of the uncertainty is given (in percentage) relative to 
the number of events in the corresponding background or signal sample. For the background, 
the size of the uncertainty relative to the total number of background events is also shown. If 
the relative uncertainties differ significantly among background samples because of statistical 
fluctuations due to the limited number of events available for the evaluation of the systematic 
uncertainties, the range from the minimum to the maximum relative uncertainty is shown. The 
dominant experimental uncertainty for the background prediction is due to the normalization 
scale factors applied to the rare EWK and V 7 samples. For the rare EWK sample, a 50% uncer¬ 
tainty is assigned as a conservative approximation of the uncertainty in the NLO-to-LO cross 
section ratio of tt 7 production, which is the dominant component in this sample. Also, for 
the rare EWK sample, we evaluate the uncertainty due to the luminosity determination H59l . 
Normalization uncertainties in the background estimates of events with misidentified photons 
or leptons are absorbed in the uncertainty of the V 7 scale factor through the uncertainty esti¬ 
mation in the A (p(£, p™' ss ) template fit described above. Subdominant systematic uncertainties 
arise from potential mismodelling of the shapes of the V 7 and proxy samples for misidentified 
photons and leptons. 

For simulation-based background estimates, differences with respect to the data in the jet en¬ 
ergy scale (JES) and resolution (JER) are considered as systematic uncertainties. The uncertainty 
associated with the JES is evaluated by varying the scale by ±1 a, where a is the half-width of 
the 68 % confidence interval around the nominal value, and recalculating the E™ ISS , Mj, and 
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Table 2: Summary of systematic uncertainties. The third column gives the uncertainty relative 
to the number of events in the corresponding background or signal sample. The fourth column 
shows, for the background terms, the uncertainty relative to the total number of background 
events in the signal region. 


Name 

Sample 

Relative uncertainty (%) 
Sample Total 

Rare background rate 

Rare EWK 

50 

19 

V 7 scale factor 

V 7 

14 

6 

Proxy sample shape 

Misidentified 7 &z£ 

20-27 

5 

Trigger and identification efficiency 

Rare EWK 

8 

3 

JES 

EWK 

0-6 

2 

V 7 shape 

V 7 

5 

2 

Integrated luminosity 

Rare EWK 

2.6 

1 

JER 

EWK 

0-2 

1 

JES 

Signal 

0-22 

— 

JER 

Signal 

0-17 

— 

Trigger and identification efficiency 

Signal 

8 

— 

Initial-state radiation 

Signal 

0-5 

— 

Integrated luminosity 

Signal 

2.6 

— 

Renormalization scale and PDF 

Signal 

4-41 

— 


Hj values in the V 7 and rare EWK samples, and similarly for the JER term. The shift in the 
expected event yield in the signal region is taken as the estimate of the systematic uncertainty. 

Table [2] also lists the systematic uncertainties considered for the signal MC samples that are 
used for the interpretation of the result of this search. The uncertainties due to the JES and 
JER are evaluated using the procedure described above. In addition, for the signal samples, 
uncertainties in the description of initial-state radiation as well as the renormalization scale 
and PDF Il60l are considered. 

6 Results 

Figure |3] shows the observed distributions of Mj, Ey ISS , photon Ej (Ey), and Hj in the e 7 and 
f /7 channels, together with the background expectations obtained as described in Section[4j The 
ratio of the observed number of events to the total background expectation is displayed in the 
lower part of each panel. The Mj distribution includes all events that satisfy the event selection 
criteria of Table [l] except for the restrictions on Mj and E™ lss . Events in the Ey ISS distribution 
must additionally satisfy Mj > 100 GeV, and events in the Ey and hlj distributions Ey ,ss > 
120 GeV and Mj > 100 GeV. The uncertainty bands shown for the background estimates are 
the statistical and systematic terms added in quadrature. The data are seen to agree with the 
SM prediction within the uncertainties. 

Figure |4] shows a compilation of event yields compared to the total background expectations. 
To enhance the sensitivity to different SUSY scenarios, the signal region is further explored 
in bins of Ey ([40,80] and >80GeV), Hj ([0,100], [100,400], and >400GeV), and E™ 1SS (Low, 
Mid, and High, corresponding to [120,200], [200,300], and >300GeV, respectively). The data 
are found to be consistent with the background prediction in all regions. Thus no significant 
excess of events beyond the SM expectation is observed. 
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Figure 3: Distributions of (a, b) Mj, (c, d) E™ ISS , (e, f) Ej, and (g, h) Hj compared with the 
stacked background expectations for the (a, c, e, g) e 7 and (b, d, f, h) uj channels. See text 
for details of the event selections applied to these distributions. The rightmost bin of each plot 
shows the overflow, with contents that are not normalized by the bin width. Expected signal 
distributions from a GMSB model for a representative mass point (Mg, = (915,405) GeV 
are overlaid. The lower part of each panel shows the ratio of the data to the predicted back¬ 
ground. The uncertainty bands represent the statistical and systematic terms added in quadra¬ 
ture. 
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19.7 fb' 1 (8 TeV) 
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Figure 4: Event yields in all signal region bins, compared with the combined SM background 
predictions. The numbers in brackets indicate the range of the bins in the signal region vari¬ 
ables ISS , Hj, and Ej. For the E™ ISS bins. Low, Mid, and High correspond to the intervals 
[120,200], [200,300], and >300 GeV, respectively. 

7 Interpretation 

The results of the search are interpreted in terms of cross section upper limits on a GMSB model 
and two distinct simplified models. For each parameter point of the three models, a large 
number of hard-scattering simulation events are generated. These events are processed with 
a detailed fast simulation of the CMS detector response I16TI . A large number of minimum- 
bias interactions are superimposed on the hard-scattering process in order to reproduce the 
pileup profile observed in data. The event selection applied to the simulated signal events is 
identical to that applied to data, including the trigger requirements. The resulting event yields 
are weighted by correction factors to account for selection efficiency differences between data 
and simulation. 

For each mass point of the signal models, a 95% confidence level (CL) cross section upper limit 
is obtained utilizing the "LHC-style" CL S prescription l62ti64l . which calculates frequentist CL S 
limits using the one-sided profiled likelihood as a test statistic. The SM background prediction, 
signal expectation, and observed number of events in each signal-region bin of the ej and ]i 7 
channels as shown in Fig.[4]are combined into one statistical interpretation, turning the analysis 
into a multichannel counting experiment. 

7.1 Interpretation in a GMSB model 

A GMSB model with wino co-NLSPs Il32l . which contains both electroweak and strong pro¬ 
duction as the primary SUSY production process, is examined. All SUSY particles except for 
the gluino and winos are considered in the limit of very large mass values such that they do 
not participate in the interactions. In this limit, the lightest chargino and neutralino become 
purely wino-like. There is no restriction on the decays, but the gluino always undergoes at 
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7 Interpretation 


least a three-body decay and the charged wino decays to a W boson and the gravitino. The 
neutral wino decays to a gravitino and either a photon or a Z boson, with branching fractions 
dictated by the weak mixing angle and the wino mass. In the generated scans, the gluino mass 
(Mg) ranges from 715 to 1415 GeV in 50 GeV steps, and the wino mass (M^) from 205 GeV to 
[Mg — 10 GeV], also in 50 GeV steps. 

The SUSY particle spectra and branching fractions are determined using the SUSPECT 2.41 1651 
and SDECAY 1.3 If 66 l programs. PYTHIA6.4 is employed for the SUSY particle generation, de¬ 
cays, and the subsequent parton showers. The cross section for each mass point is determined 
to NLO accuracy in quantum chromodynamics using the PR0SPIN0 2.1 H67l program. This 
cross section result, along with its uncertainty, is used to derive 95% CL exclusion limits on the 
SUSY particle masses. 

Figure [5] shows the observed 95% CL cross section upper limits with the exclusion contours 
overlaid. In the figure, the dashed curves are the median and ±1 standard deviation (u) ex¬ 
pected exclusion contours assuming the nominal cross sections for the signal model. The solid 
curves represent the observed exclusion with the signal cross sections at the nominal and Her 
values. All mass points on the bottom left side of the contours are excluded. Note that the ap¬ 
proximate one standard deviation discrepancy between the expected and observed exclusion 
contours toward larger values of M vv is due to an upward fluctuation observed in a single bin 
of the signal region with Ej > 80 GeV, E“ 1SS > 300 GeV, and an intermediate bin in Hj from 
100 to 400 GeV. 

In this inclusive GMSB model, electroweak production will always take place when the wino 
is light, independent of the gluino mass. Thus the exclusion curve becomes horizontal for 
M xn ~ 360 GeV. Note that for this and for all other instances in this paper where a numerical 
result is quoted for a mass limit, the result is based on the theoretical prediction for the cross 
section minus its la uncertainty. The expected distributions for signal in Fig. [3] correspond to 
the GMSB model for the mass point (Ms, M v y) = (915,405) GeV. This mass point has com¬ 
peting contributions from strong and electroweak SUSY production and exhibits a non-trivial 
behavior in Hj as can be seen in Fig. |3jg) and (h). 

7.2 Interpretation in simplified models 

In a simplified model H68H70I . a limited set of hypothetical particles and decay chains are in¬ 
troduced to describe a given topological signature such as the Ay final state studied in this 
analysis. The production and decay amplitudes of these particles are parameterized in terms 
of the particle masses. 

The two simplified models considered are denoted the TChiWg and T5Wg models. The TChiWg 
model is initiated by the direct production of hypothetical particles x^ and x°, whose decays 
are restricted to W ± G and 7 G, respectively. The gravitino G is nearly massless as in GMSB 
models. Thus, this model can be identified with electroweak production in the GMSB wino 
co-NLSP model, depicted by the diagram in Fig. [lja), differing only in the decay branching 
fractions. The particles x^ an d X° are therefore identified with gauginos in the remainder of 
this paper. A mass range of 100 < Mx — 800 GeV is considered, where My is the degener¬ 
ate mass of the gauginos. The generation of events for the T5Wg model, corresponding to 
the diagram in Fig. |TJb), starts with the pair production of gluinos. The two gluinos undergo 
three-body decays g -A qqy* and g —> qqy°, followed by the decays of the XT and as dis¬ 
cussed above. The T5Wg samples are generated in a mass region 700 < Mg < 1400 GeV and 
25 GeV < My < [Mg — 25 GeV]. No other non-SM particle is involved in either model. 
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Figure 5: Observed and expected 95% confidence level upper limits on the cross section and 
corresponding exclusion limits for the GMSB model. 

Events for both models are generated with the MadGraph 5 1.3 program, with up to two final- 
state partons in addition to the hard interaction. The events are then interfaced to PYTHIA 6.4, 
which is used to describe the SUSY decay chains and parton showers. The neutralino-chargino 
and the gluino pair production cross sections are calculated to NLO and NLO+NLL (next-to- 
leading logarithm) accuracy II60TI . respectively, and used to derive 95% exclusion limits. 

Figure [6ja) shows the computed 95% CL cross section upper limit for the TChiWg model as 
a function of My, together with the theoretical cross section. Assuming a 100% branching 
fraction for x° -A yG, the mass range 100 < M x < 540 GeV is excluded, where the lower limit 
corresponds to the lowest My included in the scan. Assuming a more physically motivated 
branching fraction through a rescaling of the theoretical cross section by the weak mixing angle, 
the exclusion range is 100 < My < 340 GeV. The latter result is similar to the limit M v y < 
360 GeV obtained from the GMSB model with a wino-like NLSP. 

The production cross section of the T5Wg model is determined solely by Mg. Nevertheless, the 
Mg — My mass difference affects the Hj and Ej' ss spectra, resulting in nontrivial exclusion-limit 
contours in the My-Mg plane. The 95% CL cross section upper limits and exclusion contours 
for the T5Wg model are shown in Fig.[6jL). For My > 200 GeV, pair production of gluinos is ex¬ 
cluded for gluino masses below 1 TeV. For 500 < My < 700 GeV, gluinos below approximately 
1.1 TeV are excluded. 
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8 Summary 

This paper presents a search for the anomalous production of events with a photon, an elec¬ 
tron or muon, and large missing transverse momentum produced in proton-proton collisions 
at i/s = 8 TeV. The data are examined in bins of the photon transverse energy, the magnitude 
of the missing transverse momentum, and Hj, the scalar sum of jet energies. The standard 
model background is evaluated primarily using control samples in the data, with simulation 
used to evaluate backgrounds from electroweak processes. No excess of events above the stan¬ 
dard model expectation is observed. The results of the search are interpreted as 95% confidence 
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Figure 6: Observed and expected 95% confidence level upper limits on the (a) TChiWg and (b) 
T5Wg simplified models. For the T5Wg model, the 95% confidence level exclusion contour is 
also shown. 


level upper limits on the production of new-physics events in the context of a gauge-mediated 
supersymmetry breaking (GMSB) model. The GMSB model is excluded for wino masses be¬ 
low 360 GeV. Results are also interpreted in the context of two simplified models inspired by 
GMSB, denoted TChiWg and T5Wg. The TChiWg model is excluded for gaugino masses be¬ 
low 540 GeV or, if the cross sections are scaled by the weak mixing angle, below 340 GeV. The 
T5Wg model with gaugino mass above 200 GeV is excluded for gluino masses below 1 TeV. 
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